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Anisotropy analysis of EAS data in the knee region 
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Abstract. Based on MAKET-ANI EAS data the distribu- 
tions of equatorial coordinates of EAS core directions are 
obtained in the knee region. Anisotropy of primary cosmic 
rays is displayed only by declination equatorial coordinates 
(Sa — 20'-* ± S'') at primary energies more than 5-10 PeV. 
The fraction of anisotropic component turns out ~ 10% in 
the knee region. 



1 Introduction 

The behavior of EAS size spectra in the knee region points 
out a possibility of a multi-component nature of the primary 
nuclear flux (Linsley, 1983; Gaisser and Stanev, 1993; Bier- 
mann, 1993; Stanev et al., 1993; Erlykin and Wolfendale, 
1997; Ter-Antonyan and Haroyan, 2000). As follows from 
(Stanev et al., 1993; Ter-Antonyan and Haroyan, 2000; Ter- 
Antonyan and Biermann, 2001), the best description of EAS 
size spectra in the knee region may be achieved providing at 
least two components in the primary cosmic ray flux. Confir- 
mation of the multi-component origin of cosmic rays can be 
obtain by investigations of distribution of EAS arrival direc- 
tions. Especially, it is interesting to measure of the anisotropy 
in the vicinity of the knee region (10^^ — 10^^ eV) where 
the accuracies of EAS experiments in last years have signifi- 
cantly increased. 

In this paper, based on MAKET-ANI (Chilingaryan et al., 
1999) EAS data 10*^ events with > 10^) the equato- 
rial coordinate distribution of EAS arrival directions are in- 
vestigated. Violations of isotropy 5 — 10%) are obtained 
at energies E > 10 PeV. 



2 EAS anisotropy 

Principal complications of the investigation of the EAS ani- 
sotropy are: 



1 . EAS attenuation in the atmosphere (the intensity of shower 
size F{Ne, 0) strongly depends on a zenith angle of in- 
cidence (6) at the given observation level); 

2. dependence of the EAS size on the energy (E) and pri- 
mary nucleus (A); 

3. interruptions of exposition local time, which are inevitable 
in EAS experiments. 

Factors 1-3 do not allow in principle to apply the traditional 
(on-off) methods of measuring the cosmic -ray anisotropy and 
here, a simulation method is applied taking into account above 
factors in the frameworks of QGSJET interaction model 
(Kalmykov and Ostapchenko, 1993) and predictions of the 
2-component origin of cosmic rays (Biermann, 1993). The 
test of this models by modern EAS size spectra at different 
zenith angles and different observation levels one can find in 
papers (Ter-Antonyan and Haroyan, 2000; Ter-Antonyan and 
Biermann, 2001). 

Assume that the primary nuclear flux I{E) consists of isotropic 
Io{E) and anisotropic components (Berezinsky et al., 

1984) 



IiE) = Io{E) + Q{E,aa,Sa) 



(1) 
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where aa , Sa are equatorial coordinates of anisotropic direc- 
tion. Let also the factor of anisotropy e = '^{E)/ I{E) is 
small (e <C 1) and does not distort the EAS zenith angular 
distribution due to the Earth rotation and independence of 
EAS attenuation on the azimuthal angle of incidence ((p). 

Then, providing each detected EAS event from I{E) flux 
with horizontal coordinates 9, (p and local time T by ade- 
quate simulated event but from isotropic Io{E) flux with 
corresponding coordinates 6* , ip* and the same T one can 
compare (after transformation 9, ip, T and 9* ,(p* ,T ^ a,d 
and a* ,6*) the obtained equatorial coordinate distributions 
of real EAS data and isotropic simulated data. Evidently, if 
the anisotropic part of intensity ^{E, aa, 5a) = then the 
distribution obtained from the experiment must overlap with 
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Fig. 1. Equatorial coordinate distributions of all showers at different primary energies (E). Symbols are the Maket-ANI data. Solid lines are 
the expected distributions obtained by simulations of the isotropic primary component (off-component) taking into account the local time of 
each detected EAS event. Dashed lines are the isotropic components at the uniform local time distribution of each EAS detected event. 



the simulated one. Any considerable discrepancies (out of 
statistical errors) will point out to the existence of anisotropy. 
The presented method completely solves the problems 1 and 
3, however, the same normalization of real and simulated dis- 
tributions (due to a flux aa, So) is unknown) does not 
give a direct possibility to determine the absolute intensity of 
the anisotropic component. 

Taking the above into account, each detected EAS event is 
characterized by a 9-dimensional vector D{E, 9, ip, T{y, m, 
d, h, u, s)), where E = E{N^, 9) is the evaluation of a pri- 
mary energy on the basis of the measured EAS size Ng and 
zenith angle 9, T{y, m, d, h, u, e) is a local time vector with 
components y -year, m -month, d -day, h -hour, u -minute, 
e -second of a detected event. Here we used the vector defi- 
nition only for briefness and convenience. 

The estimation of the primary energy E by the measured 
EAS size Ng is performed by inverse interpolation of func- 
tion = Nf.{E, A, 9) at average primary nucleus A{E) = 
exp(ln A). A values are determined by 2-component primary 
spectra (Biermann, 1993). A tabulated function Ne{E, A, 9) 
at given E, A, 9 parameters of a primary nucleus was pre- 



liminary calculated by means of EAS simulation using the 
CORSIKA code (Heck et al., 1998) at QGSJET interaction 
model (Kalmykov and Ostapchenko, 1993). Some details of 
this simulation one can find in the paper (Ter-Antonyan and 
Haroyan, 2000). Statistical eiTors did not exceed 3 — 5% and 
the high altitude location (700g/cm^) of the ANI experiment 
(ChiHngaryan et al., 1999) aflowed to obtain l\E < 20%. 

Because the exposition time of the MAKET-ANI array 
is not continuous and the EAS zenith angular distribution 
does not adequate to the angular distribution of primary nu- 
clei (factors 1,3 above), we created i — 1, . . . , m simulated 
9-dimensional vectors Ri [E* , 9* , ip* , T) for each detected 
vector D which differ from D only by simulated primary 
energy (E*), zenith (9*) and azimuth (ip*) angles. 
The values of angles 9* and ip* are simulated according to 
distributions: dF/dcos9 — cos 9 and d^ / d(p = 1/211 
spectively. The energy {E*) of simulated events is obtained 
from 2-component primary energy spectra (Biermann, 1993; 
Ter-Antonyan and Biermann, 2001) at additional condition 
of detected events N* {E* ,A*,9*)> 10^. 

Transformation of D and Ri vectors to simple 3 -dimensional 
vectors d{E,a,S) and ri{E*,a*,S*) with a -right ascen- 
sion and S -declination equatorial coordinates are performed 
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Fig. 2. Detected normalized all-particle primary energy spectra at different equatorial declination 5-coordinates. Black symbols correspond 
to AM EAS data. Open circles are the expected energy spectra of an isotropic component (off-component) according to (Biermann, 1993) 
model. 



by astronomical soft (Super- Kamiokande, 1999). 
3 Results 

The results of investigation of the EAS anisotropy in the knee 
region by AMI EAS data are shown in Fig. 1,2. Equatorial 
coordinate distributions dN/da and dN/dS of all showers 
(Ntot = 10^ events at A^e > 10^) at six primary energy (E) 
intervals are given in Fig. 1. Symbols are the Maket-ANI 
data dN/da and dN/dS in terms of (100 /Ntot)- SoHd Hues 
are the expected adequate distributions dN/ da* and dN/dS* 
obtained by simulations of isotropic primary component (off- 
component). Local times T of simulated events and detected 
EAS events are equal. The total number of events is the same 
(Ntot = ^tot)- Dashed lines reflect the expected behavior of 
the isotropic component at an uniform local time distribution 
of each EAS detected event. 

It is seen, that the distribution of a (5-coordinate at all primary 
energies is practically independent on a local time distribu- 
tion. Moreover, a good agreement of expected and detected 
distributions for equatorial a-coordinate is observed in the 



energy range -E > 0.7 PeV. 

However, 3ie discrepancies of detected and expected distri- 
butions of the equatorial (5-coordinate increase with the en- 
ergy. Strong disagreement is observed at -E > 10 PeV and 
S ~ 20'^. It is necessary to note, that the real disagreement 
is a little larger than it is seen in Fig. 1 because EAS and 
simulated data have the same normalizations. Further inves- 
tigations of the ^-coordinate distribution are given in Fig. 2 
where detected normalized all-particle primary energy spec- 
tra at different declinations (6) are presented. Black symbols 
are AMI data. The open circles are expected energy spectra 
of isotropical simulated components (off - component) with 
50 times higher statistics (N^^^ = 5 • lO"^). 
Spectral shape in E < 1 — 2 PeV range for all declination 
intervals is determined by the lower limit of the detected 
EAS size (Ng > 10^). However, the observed fine struc- 
ture of all-particle spectra at £ > 2 PeV and declinations 
S < 30" and 50 < (5 < 60 can hardly be explained by triv- 
ial statistical fluctuations. All data in Fig. 1,2 are obtained in 
the framework of QGS JET interaction model (Kalmykov and 
Ostapchenko, 1993) and the 2-component origin of cosmic 
rays (Biermann, 1993; Ter-Antonyan and Biermann, 2001). 
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4 Conclusion 

From above analysis it follows that there is a partial anisotropic 
flux with a declination coordinate equal to 5a = 20° ± 3° at 
primary energies ii^ > 10 PeV. The distribution of equatorial 
a-coordinates agrees with the hypothesis of the isotropic dis- 
tributions of cosmic ray. 

It is interesting to note, that the nearest neutron stars PSR0950 
and PSR1133 disposed at I ~ 30 — 100 ps distance from 
the Earth have declination coordinates 5 = 30.6° and S = 
27.45° respectively. These values are close to the results ob- 
tained here (da ~ 20°). Are these neutron stars a reason of 
anisotropic flux or not can be elucidated only with signifi- 
cantly large EAS statistics. 
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